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ABSTRACT 

We have imaged CO(J=7-^6) and Ci( 3 P2^ 3 Pi) emission in the host galaxy of the z=6.42 quasar 
SDSS J114816.64+525150.3 (hereafter: J1148+5251) through observations with the Plateau de Bure 
Interferometer. The region showing CO(J=7-^6) emission is spatially resolved, and its size of 5kpc 
is in good agreement with earlier CO(J=3^2) observations. In combination with a revised model of 
the collisional line excitation in this source, this indicates that the highly excited molecular gas traced 
by the CO J— 7-^6 line is subthermally excited (showing only 58±8% of the CO J=3^2 luminosity), 
but not more centrally concentrated. We also detect C i( 3 P2— > 3 Pi) emission in the host galaxy of 
J1148+5251, but the line is too faint to enable a reliable size measurement. From the Ci( 3 P2^ 3 -Pi) 
line flux, we derive a total atomic carbon mass of Mci=l.lxlO 7 M , which corresponds to ~5xl0 -4 
times the total molecular gas mass. We also searched for H20( Jif a if c =2 12 ^loi) emission, and ob- 
tained a sensitive line luminosity limit of L' n o <4.4xl0 9 K kms" 1 pc 2 , i.e., <15% of the CO(J=3^2) 
luminosity. The warm, highly excited molecular gas, atomic gas and dust in this quasar host at the 
end of cosmic reionization maintain an intense starburst that reaches surface densities as high as 
predicted by (dust opacity) Eddington limited star formation over kiloparsec scales. 
Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Detailed studies of individual quasars at the highest 
redshifts are vital to shed light on how today's most 
massive galaxies formed and to investigate whether or 
not the tight correlation between black hole mass and 
stellar bulge mass seen today was already in place dur- 
ing the earliest stages of galaxy formation. A particu- 
larly insightful example is the very high-redshift quasar 
J1148+5251 (2=6.42), which is observed at the end of 
cosmic reionization, only ^870 Myr after the Big Bang. 8 
Since its initial disco very in the Sloan Digital Sky Sur- 
vey (SDSS; Fan et al. 2003), it has been found to be an 
optical point source even at the resolution achie yed b y 
the Hubble Space Telescope (HST; White et al. l200l . 
with ultraviolet to infr ared colors typical for a type 1 
quasar (e.g., Jiang et al. 2006). From its bright and broad 
Mgn (AA2796,2803A) emission line, a supermassive black 
hole mass of 3xl0 9 M Q has been derived (Willott et al. 
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2003). It also hosts a 4.2xlO 8 M reservoir of dust that 
exhibits a high far-infrared (FIR) continuum luminosity 
of Lfir=2.2x10 13 L q (Bertoldi et al. I2003al : Beelen et 
al. 1200(1 . 

Subsequently, Jl 148+5251 was detected in mole cular 
gas emission ( carbon monoxide, CO; Walter et al. I2003h 
Bertoldi et al. l2003bfl , revealing a massive molecular gas 
reservoir of 2.4xlO lo M0 that could fuel star formation 
and feed the active galactic nucleus (AGN) in this sys- 
tem. High resolution follow-up observations revealed 
that this gas is centered on t he AG N, but distributed 
on scales of 5 kpc (Walter et al. [2004) . This large molec- 
ular reservoir harbors a compact 1.5 kpc size region that 
emits bright emission from the [C n]( 3 -P3/2— inter- 
stellar medium (ISM) cooling line, which is presumably 
due to active star formation at an enormous star forma- 
tion rate (SF R) surfa ce density of ~1000M© yr^ 1 kpc -2 
(Walter et al. I2009af ). The dense, star- forming molecu- 
lar gas component as traced by hydrogen cyanide (HCN) 
is comparatively faint, given th e sou rce's high SFR and 
FIR luminosity (Riechers et al. 120071 ). 

In addition, it was found that J1148+5251 is radio- 
quiet and follows the radio-FIR correlation for star- 
forming galaxies, providing additional evidence that the 
FIR dust emission is dominantly hea ted b y young stars 
rathe r than the AGN (Carilli et al. 120041: Wang et al. 
2008). To further investigate on which scales the star 
formation takes place, we here aimed at resolving the 
neutral atomic and the warm, highly excited molecular 
gas components in the host galaxy of this unique quasar. 



2. OBSERVATIONS 
We observed the CO(J=7— >6) (f r cst 
806.6518GHz), Ci( 3 F2^ 3 -Pi) (809.3435 GHz), 
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RIGHT ASCENSION (J2000) 



Fig. 1.— Velocity-integrated PdBI map of the CO(J=7— >6) 
line emission toward J1148+5251 (contours) over the central 
41 4 km s -1 , overlayed on the VLA CO(J=3^2) map (Walter et 
al. 120041 : color scale). At a resolution of 0.86"x0.61", the emission 
is resolved. The cr oss indicates the position of the optical quasar 
(White et al . 12003) . Contours are shown at (-3, -2, 2, 3, 4, 5, 6)xcr 
(lo- = 0.125 mjy beam" 1 ). 

H 2 0(Jk q k c =2i2^1oi) (1669.9048GHz) emission lines 
toward J1148+5251 using the IRAM Plateau de Bure 
Interferometer (PdBI). At z=6.419, these lines are 
redshifted to 108.7278, 109.0907, and 225.0849 GHz (2.8 
and 1.3 mm). First observations were carried out with 
the PdBI's previous generation receivers in the old 6D 
configuration in May, June and Novem ber 20 03 (May 
observations published by Bertoldi et al. |2003bl ), and in 
April-June and September-November 2004 (23 tracks 
total). Further observations were carried out with the 
new generation receivers in the new 6A configuration 
(longest baseline: 760 m) in February and March 2007 
(7 tracks total). All tracks were taken under good 3 mm 
observing conditions. Of the tracks taken with dual 
frequency setup in 2003/2004, 12 were taken under 
conditions sufficient for 1 mm observing. The nearby 
source 1150+497 (distance to J1148+5251: 3.4°) was 
observed every 22.5 minutes for pointing, secondary 
amplitude and phase calibration. For primary flux 
calibration, several nearby standard calibrators were 
observed during all runs, leading to a calibration that 
is accurate within 10-15% (3mm) and 15-20% (1mm), 
respectively. 

Observations with the previous generation receivers 
were set up using a total bandwidth of 580 MHz 
(single polarization, dual frequency; corresponding to 
-1600 kms -1 at 2.8mm and -770 kms -1 at 1.3mm). 
This setup requires two separate tunings at 3 mm to 
observe the CO(J=7->6) (14tracks) and Ci( 3 P 2 ^ 3 Pi) 
(9 tracks) emission lines (1mm receivers were tuned to 
the redshifted H 2 J K a K a — ^12^^01 frequency). Ob- 
servations with the new generation receivers were set 
up using a total bandwidth of 1 GHz (dual polarization; 
corresponding to —2800 kms' 1 at 2.8mm). This setup 
allows to observe the CO(J=7-*6) and Ci( 3 P 2 ^ 3 Pi) 
lines simultaneously (selecting a tuning frequency of 
108.894GHz centered between both lines). All observa- 
tions were taken with sufficient bandwidth to cover the 
underlying continuum simultaneously. 
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Fig. 2. — PdBI spectrum of the CO(J=7— »6) and underlying 
2.8 mm continuum emission toward J1148+5251 at a resolution of 
28 kms -1 (10 MHz). The velocity scale is relative to the redshifted 
CO(J=7— >6) frequency at 2=6.419. The rms per velocity bin is 
0.33 mjy. The solid line shows a Gaussian fit to the data. 

For data reduction and analysis, the IRAM GILDAS 
package was used. All data were mapped using 'natural' 
weighting unless mentioned otherwise. The CO( J=7^6) 
data result in a final rms of 0.33 mjy beam -1 per 
28 kms -1 channel (resolution: 1.38"xl.l7"). To opti- 
mize the spatial resolution, these data were also mapped 
using 'uniform' weighting, leading to a synthesized clean 
beam size of 0.86" x 0.61" (4.8 kpc x 3.4 kpc) and an 
rms of 0.12 mjy beam -1 over 414 kms -1 (150 MHz). To 
maximize the signal-to-noise ratio, the Ci( 3 P 2 ^ 3 Pi) 
data were tapered to a resolution of 2.78" x 2.42", lead- 
ing to an rms of 0.17mJy beam -1 over 316 kms -1 
(115 MHz). The 2.8mm continuum data (combined from 
all CO J=7— >6 and C i 3 P 2 — > 3 Pi observations) result in a 
1.20"x0.99" beam and an rms of 0.046 mjy beam -1 over 
563.75MHz. The H 2 0( J KaKc =2 12 ^l 01 ) data result in 
an rms of 0.76 mjy beam -1 over 303 kms -1 (227.5MHz), 
and the (double sideband) 1.3 mm continuum data result 
in an rms of 0.34 mjy beam -1 over 2x563.75 MHz. 

3. RESULTS 

In Figure[I] the velocity-integrated CO(J=7— >6) emis- 
sion is shown at a linear resolution of —4 kpc (0.7"). 
From elliptical Gaussian fitting to the Her detection of 
the source in the u — v plane, it is found that the emission 
is spatially resolved on a scale of 0.9"±0.16" (5.0 kpc) 
in the north-south direction, and marginally resolved 
(0.54"±0.11"; 3.0 kpc) in east-west direction. The peak, 
extent 9 and orientation of the emission are in remark- 
able agreement with the lower-excitation CO( J =3— >2 ) 
line emission (color scale in Fig. Q] Walter et al. 120041 ). 
and the peak is coincident with the positi on of the op- 
tical quasar (cross in Fig. [Q White et al. [20051) . This 
demonstrates that the relative astrometry between the 

9 Accounting for beam convolution. 
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Fig. 3. — Map of the 2.8 mm continuum emission toward 
J1148+5251. At 1.20" xO. 99" resolution, the emission appears 
unresolved. Contours are shown at (-3, -2, 2, 3, 4, 5)xcr (1<t = 
0.046 mjybeam" 1 ). The cross indicates the same position as in 

Fig. m 

CO(J=7^6) PdBI and CO(J=3^2) VLA observations 
is accurate within <0.1". It also indicates that, even 
though the CO(J=7^6) line is only subthermally ex- 
cited (see discussion below), the emission is apparently 
not more centrally concentrated than the CO(J=3^2) 
emission. 

In Figure[2j the spectrum of the CO( J=7— >6) emission 
is shown. The line is detected at a peak flux of S u = 
1.99 ± 0.23 mJy and a width of 297±35 kms" 1 , leading 
to an integrated line flux of 0.63±0.06 Jy kms -1 , and a 
line luminosity of L' co(7 _ 6) =1.7+0.2 x 10 10 K kms^pc 2 

(see Table []} . The sensitivity of the observations is suffi- 
cient to detect, for the first time, the underlying 2.8 mm 
continuum emission (0.28±0.07mJy). 

Figure [3] shows a map of the 2.8 mm continuum emis- 
sion integrated over all line- free channels of the combined 
CO and Ci observations at ~6kpc (1.1") linear resolu- 
tion. The continuum emission is detected at a flux level 
consistent with the simultaneous line/continuum fit to 
the CO(J=7^6) spectrum within the errors. 

In Figure [4j a velocity-integrated map of the 
Ci( 3 P2^ 3 Pi) emission over 316 kms -1 (i.e., approxi- 
mately the CO line FWHM) is shown at a linear res- 
olution of ~15kpc (2.6"). Emission is detected at 
0.98±0.17mJy peak flux, corresponding to a line peak 
flux of 0.7 mJy when accounting for the underlying con- 
tinuum. This corresponds to an integrated C i( 3 P2— > 3 Pi) 
line flux of 0.22±0.05 Jy kms -1 , and a line luminosity 
of L' CI , 2 ^=6.0±1.3xl0 9 K kms -1 pc 2 . At higher spa- 
tial resolution, the line is still marginally detected, but 
the peak flux decreases. This may indicate that the 
C i( 3 p2 — >3 Pt) emission is resolved on similar scales as the 
CO emission, and thus emerges from the same (molecu- 
lar) gas phase on global (galactic) scales. However, this 
conclusion remains tentative at the present signal-to- 
noise ratio. 

No evidence for H20( Jn a K c =2i2^1oi) emission is 
found in the 1.3 mm data. Assuming that the line has 
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Fig. 4. — Velocity-integrated map of the C l( 3 P2 _ * 3 -Pi) line emis- 
sion over 316 kms -1 toward J1148+5251. At 2.78"x2.42" resolu- 
tion, the emission appears unresolved. Contours are shown at (-3, 
-2, 2, 3, 4, 5)xct {la = 0.167mjy beam" 1 ). The cross indicates 
the same position as in Fig. [T] 

the same width as CO(J=7— >6), we derive a 3(7 upper 
limit of <0.69 Jy kms -1 for the integrated line flux, and 
a line luminosity limit of 

^ 2O <4.4xl0 9 Kkms- 1 pc 2 (i.e., £^0/^0(3-2) < 
0.15). We clearly detect the underlying 225 GHz con- 
tinuum emission at 3.9±0.8mJy. 

4. DISCUSSION 

4.1. Properties of the Molecular and Atomic ISM in 
J1148+5251 

4.1.1. Resolved CO(J=7—y6) Emission 

We successfully resolved emission from the high- J 
CO( J=7^6) line toward the z=6.42 quasar J1148+5251. 
The structure and size of the gas reservoir are consis- 
tent with what was derived from previous CO(J=3^2) 
line mapping (Walter et al. 12004ft . The line width de- 
rived from our high signal-to-noise CO(J=7^6) spec- 
trum is consistent with that of previous observations of 
the CO(J=6-^5) and CO(J=7->6) lines (Bertoldi et al. 
2003b). Although the molecular gas emission appears 
to be more extended than that from the neutral ISM 
as traced by the 158 /im [C n]( 3 P 3 /2^ 3 Pi/2) line, the 
linewidths and centro ids are consis tent to high precision 
(Maiolino et al. 120051 : Walter et al. l2009an . 

4.1.2. Atomic Carbon 

We also detected the upper fine structure line of neu- 
tral carbon, Ci( 3 P 2 ^ 3 Pi), toward J1148+5251. This 
represents the highest z detection of atomic carbon (or, 
indeed, any neutral atomic medium) to date, and the 
first C i detection in an unlensed high-rcdshift galaxy. 
As the lower fine structure line (which is redshifted to 
66.3 GHz) cannot be observed with current facilities, it 
is not possible to derive an excitation temperature for 
Ci. However, for T ex > 20 K, the de rived C i mass only 
weakly depends on T ox (Weifi et al. 2005a). Assuming 
an excitation temperature of T ex =36K (the average of 
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TABLE 1 

Line fluxes and luminosities in SDSS ,11148+5251. 



S„dv V Ref. 

[Jykms- 1 ] [10 9 K kms" 1 pc 2 ] 



CO(J= 


=1-0) 


<0.11 


<143 


1 


CO(J= 


=3—2) 


0.20 ± 0.02 


29.9 ± 2.7 


2 


CO(J= 


=6—5) 


0.67 ± 0.08 


25.0 ± 3.0 


1 


CO(J= 


=7—6) 


0.63 ± 0.06 


17.2 ± 1.8 


1,3 


HCN(J= 


=2-1) 


<0.006 


<3.3 


4 


HCO+(J= 


=2-1) 


<0.018 


<10 


5 


H 2 0(JK a jf c =2i2 


-loi) 


<0.69 


<4.4 


3 


Cl( 3 P 2 


- 3 Pi) 


0.22 ± 0.05 


6.0 ± 1.3 


3 


[CH]( 3 P 3 /2^ 




3.9 ± 0.3 


19.0 ± 1.6 


6,7 


[NH]( 3 P r 


- 3 Po) 


<0.47 


<4.0 


8 



4.1.3. Limits on Water Emission 

We have also searched for H 2 0( J#: i<: e =2i2— ►loi) emis- 
sion toward J1148+5251, but only obtained a sensi- 
tive upper limit. Radiative transfer calculations show 
that, assuming the known gas and dust properties 
of J1148+5251, H 2 0( JK a A" c =2i2^1oi) emission is ex- 
pected to be optically thick over a substantial range of 
the permitted parameter space. Thus, the line luminos- 
ity limit can be directly translated to a surface filling 
factor limit (relative to CO J=3^2) of <15% (assuming 
both lines are in thermal equilibrium). 

Sensitive searches for water emission from dense molec- 
ular gas were perfor med for two othe r z>3 galaxies 
(Riechers et al. I2006at Wagg et al. 2006) , reaching com- 
parable depths. A search for H 2 0( Jk q a: c =3i3— 2 2 o) 
emission (t'rcst^lSS.SlOGHz) toward the z=3.20 quasar 
MG 0751+2716 provided a surface filling factor limit 
of <6% relative to CO(J=4 — >3) (a ssuming the same 
conditions; Riechers et al. 2006a). A search for 
H20(JK a K c =lio^loi) emission (i/ rest =556.936GHz) to- 
ward the 2=3.91 quasar APM 08279+5255 provided 
a surface filling factor limit of <12 % relat ive to 
CO( J=4- ± 3UWagg et al. [20061 WeiB et al. l2007at Riech- 
ers et al. 120091) . 

Note that the H20(JK- if =3i3— ^20) transition was 
detected toward the nearby ULIRG Arp 220, at about 
one third of its H CN(J=1— 0) line luminosity (Cer- 
nicharo et al. [2006) . This suggests a small surface filling 
factor even within dense molecular cloud cores as traced 
by HCN. This is, again, assuming optically thick emis- 
sion in thermal equilibrium from both lines. Note that 
the H2O line emission in Arp 220 (and the high-z sources) 
is likely maser-enhanced, which would even predict sub- 
stantially smaller su rface fi lling factors (Cernicharo et al. 
120061 : Riechers et al. l2006aD . Although different H 2 and 
HCN transitions were observed toward J1148+5251, the 
observed H20( JK a K c — ^12— >loi) luminosity limit would 
be consistent with a H2O/HCN line luminosity ratio sim- 
ilar to Arp 220. 

4.2. Continuum Emission 

The new detections of the observed-frame 225 GHz 
(rest-frame: 1.67 THz) and 109 GHz (807 GHz) contin- 
uum emission are consisten t with the overall dust SED 
of J1148+5251 (Beelen et al. 120061) . In particular, the de- 
tection of 2.8 mm continuum emission allows us to bet- 
ter constrain both the CO(J=7— 6) and CO(J=6— 5) 
line fluxes. From the shape of the SED, we derive 
an estimated continuum flux of 0.19 mJy at the wave- 
length of the CO(J=6— 5) emission, by which we reduce 
the CO(J=6— > 5) flux relative to the original estimate 
(Bertoldi et al. I2003bl ) in the following. 

5. MODELING AND CONCLUSIONS 

5.1. CO Line Excitation 

Given the high signal-to-noise of the CO( J=7^6) ob- 
servations and the fact that we successfully detected the 
(observed-frame) 2.8 mm continuum emission (which also 
provided an updated CO(J=6— >5) line flux), we have 
calculated new Large Velocity Gradient (LVG) models 
of the line excitation in this system (see Bertoldi et al. 
2003b for the original study). In this LVG model study, 
the kinetic gas temperature and density are treated as 



References. — [1] Bertol di et a l. I l2003bl) , [2] Wal ter et al. 112003) , [3] 
this w ork, [4 ] Riechers et al. 1 120071). [5 1 Carilli et al. (120051). [6| M aiolino 
et al. d2005fl . [7] Walter et al. <2009an . [8] Walter et al. i200m) . 

Note. — CO(J=3— 2) corrected by 10% to account for flux in 
the outer linewings not covered by the observations. Line fluxes are 
continuum-corrected where applicable. 

the Cloverleaf and IRAS F10214+4724, the only high-z 
galaxies for whi ch bot h transitions have been de tected; 
Ba rvaini s et al. [19971 : WeiB et al. l200l I2005at Ao et 
al. 120081 ). we derive a Ci mass of M C i=l.lxlO 7 M . 
From the CO luminosity, a total molecular gas mass of 
M H2 =2.4xlO lo M can be derived (Walter et al. [200l . 10 
We thus find a C1/H2 mass fraction of 5xl0 -4 . Assum- 
ing a [C1I/IH 2I abundance ratio of 5xl0~ 5 (WeiB et al. 
2003, 2005a), the C 1 mass translates to a total molec- 
ular gas mass of Mh 2 =3.7x10 10 M q . This is somewhat 
higher than the value derived from the CO luminosity, 
and may indicate either a higher atomic carbon abun- 
dance, or that the CO luminosity to H2 mass conversion 
factor is higher than that assumed above [possibly due to 
different molecular abundances relative to nearby ultra- 
luminous infrared galaxies (ULIRGs)]. 

By comparing the line luminosities (in units of Lq), 
one finds that the Ci( 3 P2^ 3 Pi) line provides only ^2% 
of the cooling capacity of the [C 11] ( 3 p3/2^ 3 -Pi/2) line 
(Maiolino et al. 120051 : Walter et al. l2009af >: the latter thus 
clearly remains the strongest detected coolant of the neu- 
tral atomic ISM. This also implies Lci/iFiR — 5 x 10~ 6 . 
This ratio is comparable to those found in nearby star- 
forming and starbursting galaxies, where it is considered 
to be a measure for the strength of the non-i onizin g stel- 
lar U V rad iation field (e.g., Gerin & Philips 120001 : Bayet 
et al. [2006) . This supports the assumption that the bulk 
of the dust and gas heating in this system is powered by 
star formation. Together with the I/jctti/I/fir ratio of 
~2x 10~ 4 (Walter et al. l2009af ), the above ratios imply 
that the UV radiation field in Jl 148+5251 is relatively 
strong, and that the gas density is high [as predicted by 
photodissociation region (PDR) models] , con sistent with 
nearby ULIRGs like Arp 220 (Gerin & Philips [2000J, their 
Fig. 8). 

10 CO in high-z quasar hosts is more widespread than in nearby 
ULIRGs, but shows similar physical properties. We thus adopt 
a low ULIRG CO luminosity to H2 mass conver sion f actor of 
o=0. 8M (K kms -1 pc 2 ) -1 (Downes & Solomon 1998} rather 
than a=4-5 Mp, (K kms -1 pc 2 ) - 1 as i n near by spirals (e.g., Scov- 
ille & Sanders [19871; Solomon & Barrett HMO) . S uch low a are also 
found for 2 ~ 2.5 submm galaxies (Tacconi et al. 120081) . 
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Fig. 5. — Left: CO excitation ladder (points) and LVG model (line) fo r J1148 +5251. The upper limit for CO( J=l— >0) and the (continuum 
corrected) data point for CO(J=6— >5) arc taken from Bertoldi et al. (2003b). The CO(J=3— >2) data point (corrected for m issing fl ux in 
the line wings ) is taken from Walter et al. <2"003ft . The grey points indicate the original, uncorrected fluxes from Bertoldi et al. fl2003bT ). The 
model predicts kinetic temperatures and gas densities of Tk; n = 50 K and p ga s(H2) = 10 4 ' 2 cm - 3 . Right: Comparison of the normalized 
excitation ladders for different sources: the inner disk of the Milky Way (tr iangles and thick dashed line; Fixsen et al. I1999Q , the z=2.5 
submillimeter galaxy SMM J16359+6612 (thick solid line; Weifi et al. 20 05d ), the ne arby starburst galaxy NGC 253 (dashed line; Glisten 
et al. l2006ft . the z=4.7 quasar BR 1202-0725 (dash-dotted line; Riechers et al. l2006bn . and J1148+5251 (squares and solid line). 

free parameters. For all calculations, the H2 ortho-to- 
para ratio was fixed to 3:1, the cosmic microwave back- 
ground temperat ure w as fixed to 20.25 K (at 2=6.42), 
and the Flower (2001) CO collision rates were used. 
We adopted a CO abundance per velocity gradient 12 of 
rC OI/fdWdr) = 1 x 10- 5 pc(kms-J_r 1 (e.g., Weifi et 
al. I2005U I2007at Riechers et al. I2006bft . and a CO disk 
radius of 2.5 kpc as found from the spatially resolved 
CO(J=3->2) and CO(J=7^6) observations. The best 
solution was obtained for a spherical, single-component 
model with a CO disk filling factor of 0.16, Tk; n =50K, 
and /9 gas (H 2 ) = 10 4 - 2 cm -3 (Fig. [5]). These values are com- 
parable to th ose fo und for ot her high -z quasars (e.g., 
Riechers et al. l2006bt Weifi et al. l2007bl ). This model pre- 
dicts that the CO(J=3— »2) line emission is thermalized. 
Also, the model-predicted gas temperature is consi stent 
with that of the dust (55 ± 5 K; Beelen et al. 120061) . 

5.2. Star Formation Rate and its Surface Density 

The dust reservoir in Jl 148+5251 exhibits a FIR con- 
tinuum luminosity o f Lftr =2.2 x 10 13 Lq (Bertoldi et al. 
I2003al Beelen et al. 120061 ) . about half of which emerges 
from a compact, 0.75 kpc radius region (Walter et al. 
I2009al) . Its molecular gas, dense gas, radio continuum 
an dust properties are consistent with Lfir being dom- 
inantly powered by star formation (e.g., Beelen et al. 

1 This is due to the relative statistical weights of the symmet- 
rical (ortho) and antisymmetrical (para) eigenstates of the wave- 
function: there are three symmetrical combinations of the spins of 
both H nuclei, but there is only one antisymmetrical combination. 

12 Only the ratio of [CO] and dv/dr enters the LVG calculations. 
Thus, our solutions can account for the factor of a few variations 
in [CO] found between nearby star-forming galaxies by adjusting 
the 'internal' parameter dv/dr accordingly. Here, [CO]/(di)/dr) is 
fixed to that of the nearby starburst M82 (Weifi et al. 2005bh. 



[20061: Riechers et al. 120071 Walter et al. I2009al ). As- 
suming the AGN contribution to Lfir is small, we thus 
derive 13 an integrated SFR of 3300 M Q yr -1 , about half 
of which takes place within a 0.75 kpc radius region. 
This high SFR is consistent with the high luminosity 
of the [C n]( 3 P 3/2 ^ 3 F 1/2 ) ISM cooling line (Maiolino et 
al. l2005ft . which emerges from a region of 0.75 kp c radiu s 
within the molecular gas reservoir (Walter et al. [2009a). 
An intense, kpc-scale starburst is also needed to explain 
the fact that the molecular gas is warm and highly ex- 
cited over a large, 5 kpc size region (as shown by the new, 
spatially resolved CO J—7~ >6 data), which requires a 
heating source of similar size (given that radiation from 
a compact heating source such as the AGN likely cannot 
penetrate the dense molecular reservoir out to such large 
scales). 

Thus, the (peak) size and brightness of the resolved 
[Cii]( 3 P 3 / 2 — > 3 P 1 / 2 ) line and FIR continuum emission 
are consistent with a nuclear 'hyper'-starburst as the 
main heating source for the gas and dust, i.e., a star- 
burst that exhibits a flux of Ffir=10 13 Lq kpc -2 , pro- 
duced by an enormous SFR surfac e densi ty of Esfr = 
1000M Q yr -1 kpc -2 (Walter et al. l2009al) . Such high 
Lfir and Esfr are typically observed in the most ex- 
treme star-forming environments in the local universe. 
As a reference, the Galactic star forming cloud Orion 
exhibits a FIR luminosity of 1.2xl0 5 L Q within its cen- 
tral arcmin 2 (0.013 pc 2 ; Werner et al. 11976ft . This 
corresponds to a Lfir of ^10 13 Lq kpc -2 . Nearby 
ULIRGs exhibit a FIR luminosity of typically 3xl0 n Lq 
within their most active ~100pc size regions (Downes 

^Assuming SFR=1. 5xl0- 10 L FIR (M yr- 1 /^©) (Kennicutt 

dssi). 
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& Solomon |1998[ ). This, again, corresponds to a .Ffir 
of ~10 13 L & kpc~ 2 . Fpm ( or S S fr) as extreme as 
in J1148+5251 thus are observed in the Galaxy and 
the local universe, albeit in areas that are by 2-8 or- 
ders of magnitude smaller. Over surface areas com- 
parable to J1148+5251, extreme starbursts at z > 2 
(submillimeter galaxies; SMGs) h ave ty pical Ssfr of 
'-SOMQyr^kpc- 2 (Tacconi et al. 12006ft . which are by 
about an order of magnitude smaller. 

Intriguingly, recent high resolution molecular gas and 
radio continuum observations of the z=4.4 merging 
quasar host of BRI 1335-0417 indicate Ssfr comparable 
to those in J1148+5251 (assuming that the FIR-cmitting 
dust is n ot more extended than the m olecular gas; Riech- 
ers et al. l2008al Momjian et al. l2007f ). 

5.3. Radiation- Pressure Supported Starburst Disks 

A SFR surface density of S SF r = 1000 Mq yr" 1 kpc~ 2 
is consiste nt wit h theories of 'maximum starbursts' 
(Elmegreen 119991 ) . At such high surface densities, feed- 
back from star formation is likely to occur, which may 
(self-)regulate the starburst. If the maximum intensity 
of a radiation pressure-supported, galaxy scale starburst 
is determined by the Eddington limit for dust, a number 
of characteric limits arise in the se 'maximum starburst' 
theories (Thompson et al. 2005). Assuming a Salpeter- 
like initial mass function (IMF), a constant gas-to-dust 
ratio with radius, and that the disk is self-regulated (i.e., 
Toomre-Q~l), such an Eddington-limited starburst has 
S SFR ~ lOOOMoyr-ikpc" 2 , F ¥m ~ 10 13 L Q kp C - 2 , 
and an effective temperatur e of 8 8 K (using equations 
34-36 of Thompson et al. |2005[ ). As shown above, 
J1148+5251 has Ssfr and Ffir close to these theoretical 
limits over a 0.75 kpc radius region. 

For a flat disk as assumed here, 
-Ffir=£fir/(7t r 2 )=<7sB Trfj , where r is the disk ra- 
dius, and ctsb is the Stefan-Boltzmann constant. 14 
About half of Lfir comes from a 0.75 kpc radius region 
in J1148+5251. This corresponds to T o p ff cak =82K, which 
is also consistent with the limit. Note that, integrated 
over the whole 5 kpc region traced by CO (and full 
^fir), this simplified calculation predicts T c ff=53K, 
which is in good agreement with the kinetic temperature 
of the molecular gas as derived above, and the dust 
temperature obtained from fitting the SED of the source. 

The vertical optical depth of a disk is given by 
7V=S gas /i/2, where S gas is the gas surface density, 
and k ~ k T^ s (valid for T cS < 200 K; n ~ 2.4 x 
10 _4 cm 2 g -1 K -2 ; see Thompson et al. |2005[ ) is the 
Rosseland mean opacity for dust. The assumption of a 
constant gas-to-dust ratio with radius implies that about 
half of the gas mass is found within a 0.75 kpc radius re- 
gion in J1148+5251. Together with the T e p ff cak derived 
above, this implies ry=1.03 (i.e., moderately optically 
thick), consistent with the assumptions of the Eddington- 
limited starburst model. 

From the dynamical [Cn]( 3 P3/2— line map 
of J1148+5251, we can derive a dynamical mass of 
Miyn~l-5x 10 10 Mq within a 0.75 kpc radius region. As- 
suming that Mdyn traces the total mass within that re- 

14 Assuming a black-body, which the dust SED only approxi- 
mates. 
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Fig. 6. — Molecular gas peak surface density (Egas^) as a func- 
tion of FIR luminosity (Lfir) for nearby luminous and ultra- 
luminous infrared galaxies (squares; Wilson et al. 20(33) and 
z > 4 quasar host galaxies (crosses), including PSS J2322+1944 
(z=4.12), BRI 1335-0417 (2=4.41), BR 1202-0725 north and south 
(z=4.69), and J1148+5251 ( z=6.42; this work; Carilli et al. 12003 ; 
Riechers et al. [2008a, 2008b). Data for J2322+1944 are corrected 
for gravitational lensing. For comparison, the surface-averaged 

Ega/ 2 at half intensity radius for z > 2 SMGs is shown (circle; 
sample average of Tacconi et al. 120061) . 

gion, this corresponds to a gas mass fract ion of / gas =0.72. 
Using equation 37 of Thompson et al. (|2005f) . this sug- 
gests that the disk is optically thick out to a radius of 
^400pc. 15 This value is by almost a factor of 2 smaller 
than 750 pc. This may indicate that the potential is not 
isothermal, as assumed in the model. It may also indicate 
that the three-dimensional distribution of the gas and 
dust is more complicated than assumed here. However, 
the model reproduces a number of the observed prop- 
erties of J1148+5251 fairly well. Our observations thus 
are, indeed, consistent with a kiloparsec-scale 'maximum 
starburst', radiating close to its characteristic Eddington 
limit. This hyper-starburst is harbored by a large, 5 kpc 
size, overall warm and highly excited molecular gas reser- 
voir that hosts 2.4 x 10 10 M Q of molecular hydrogen and 
l.lxl0 7 M Q of atomic carbon. 

5.4. Dense Molecular Gas and Gas Surface Density 

There are additional lines of evidence that suggest 
that Ssfr as high as 1000 M yr -1 kpc~ 2 on kpc scales 
(as observed in J1148+5251) are not unreasonable. Re- 
cent investigations of the (in the nearby universe lin- 
ear) HCN-FIR luminosity correlation at high z (includ- 
ing Jl 148+5251) show tentative evidence for an excess in 
Z/fir relative to the observed ijjCN toward the highest 
£f tr sy stems (both with and without AGN; Carilli et 
al. (20051 : Gao et al. 120071 : Riechers et al. I2007f l. Such an 
excess can be explained by higher median gas densities 
(comparable to the critical density of HCN J=l— >0) in 

15 As Mtiy n 0C(Tp W „ M was used to determine the gas fraction, 
this value is independent of the assumed disk inclination. 
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the most FIR-luminous g alaxies (Krumholz & Thomp- 
son [2007| Narayanan et al. 2008a). Even when assuming 
a constant efficiency of turning a unit mass of molecu- 
lar gas into stars among all galaxies, this would imply 
a higher SFR per unit area for the most FIR-luminous 
systems. 

The (molecular) gas surface density in J1148+5251 

consistently reaches high peak values of up to £ g as ak ^ = 

3.5 x 10 9 M Q kpc- 2 . In Figure H the S^ ak) of 
J1148+5251, other z > 4 quasar host ga laxies (Carilli 
et al. 2002; Riechers et al. 2008a, 2008b]), and a sam- 
ple of nearby LIRGs and ULIRGs (Wilson et al. 12003) 
is shown as a function of Lfir- For comparison, a lowe r 
limit for z > 2 SM Gs is shown (Tacconi et al. [2006). 
Wilson et al. (2008) suggested that there may be a trend 

of increasing Sgas ^ with _Lfir.- This would be consis- 
tent with a higher density of gas clouds toward more 
intensely star-forming galaxies. However, their sample 
of nearby galaxies spans only about an order of magni- 
tude in Lfir- Including the more FIR-luminous high- z 
systems, this trend gets substantially weaker. In fact, 
the high-z systems appear to occupy a similar range in 

Sgas ak ' as the nearby ULIRGs, and the relation appears 
to saturate close to the typical gas surface densities of Gi- 
ant molecular clouds (GMCs). This would imply that, in 
the most FIR-luminous systems, the entire, kiloparsec- 

scale regions with the highest Sgas^ have gas surface 
densities close to those in GMCs. One limitation of this 
conclusion is that the nearby galaxies were studied at 
a factor 1.8 higher linear (physical) resolution than the 
z > 4 systems. This, however, is not sufficient to explain 
the weakening of the trend. In particular, the z=4.4 

quasar BRI 1 335-0417 with the highest S g ^ ak) (Riechers 
et al. I2008al ) was studied at only 20% lower resolution 
than the nearby sample average. Even assuming conser- 
vatively that S g as ak ^ is by ^50% higher in this source (ac- 
cording to the difference in beam area), it is still consis- 
tent with that in the nuclei of the nearby ULIRG Arp 220 
(Downes & Solomon [1998) . This would be consistent 
with a picture in which the molecular ISM in ULIRGs 
and even more FIR-luminous systems is dominant ly in a 
rather dense, relatively continuous phase (rather than in 
individual, moderately dense molecular clouds), exhibit- 
ing optically thick emission. While the measured (peak) 
gas surface densities are comparable, the gas is spread 
out over larger areas in the high-z quasar host galaxies 
relative to the (dense) nuclei of nearby ULIRGs, con- 
sistent with the higher total gas masses and SFRs, and 
higher median gas densities (integrated over the whole 
galaxies). Overall, this may suggest that a tighter cor- 
relation would be found if Xg!a| is compared to £sfr 
instead of the integrated Lfir (as a proxy of SFR). 

5.5. Evolution of the Gaseous and Stellar Components 

Assuming that the massive amounts of molecular mate- 
rial available on 5 kpc scales in J1148+5251 are converted 
into stars at an e fficienc y of 5%-10% a s in G MC cores 
(e.g., Myers et al. 119861 : Scoville et al. I1987D . the SFR 
corresponds to a gas depletion timescale of only (0.7-1.4) 
x 10 s yr. This is about one-t hird of the value found for 
BRI 1335-0417 (Riechers et al. l2008al) . suggesting a rapid 
buildup of the stellar component in both galaxies. Still, 



this gas depletion timescale is a few times larger than the 
dynamical and free-fall times of both the 0.75 kpc radius 
region traced by [C if] , and the full 5 kpc size molecular 
reservoir. In this regard, the ongoing star formation may 
still be considered slow. 

The [C n] ( 3 P 3 /2^ 3 ^i/2) dynamical mass estimate sug- 
gests that J1148+5251 hosts a total mass of ~1.5 x 
10 10 Mq within a ~0.75kpc radius region. A present 
day, massive elliptical galaxy with a black hole mass 
comparable to J1148+5251 and a velocity dispersion of 
300 kms -1 hosts a total mass of ~5-8 x 10 10 M Q within 
a central region of the same size (and a total stellar mass 
of ^2 x 10 12 Mq di stribu ted over its entire spheroid; 
e.g., Flaring & Rix 2004). Even assuming that the 
whole amount of molecular material in J1148+5251 were 
to eventually contribute to the assembly of the stellar 
spheroid within this small central region leaves a factor 
of a few difference between the central mass budgets of 
the z=6.42 quasar and a z=0 massive elliptical galaxy. 
This suggests that J1148+5251 has to accrete additional 
material by z=0 (e.g., through subsequent mergers) to 
grow a spheroid mass comparable to its (likely) present- 
day counterparts. 

On the other hand, the widths of the CO(J=7^6) and 
[Cii]( 3 P 3 / 2 — >3 -Pi/2) lines are the same within the uncer- 
tainties, while the sizes of the CO and [C njemitting re- 
gions appear to be different by a factor of ~3. Thus, [CII] 
traces a smaller dynamical mass than CO, corresponding 
to a shallower potential well. This difference may indi- 
cate that the gas in Jl 148+5251 has not yet fully coa- 
lesced, possibly due to an ongoing merger. This would be 
consistent with the high star formation efficienc y foun d 
in the central region of this galaxy (Walter et al.|2009a), 
and the compa ratively high fraction of dense gas (Riech- 
ers et al. 120071) . 

We conclude that the highly excited, several kpc scale 
size gas reservoirs in dust- and gas-rich high-z quasar 
host galaxies are dominantly heated by the large-scale 
starbursts that they maintain, consistent with cosmolog- 
ical sim ulation s of merger-driv en z ~ 6 quasar formation 
(Li et al. l2008t Narayanan et al. l2008bl ). These starbursts 
reach surface densities as predicted by Eddington- limited 
star formation over kpc scales, accompanied (and prob- 
ably supported) by ongoing major, 'wet' merger activity 
in some cases. 
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